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The Hyperspherical Harmonics basis, without a previous symmetrization step, is used to calculate 
binding energies of the nuclear A = 6 systems using a version of the Volkov potential acting only 
on s-wave. The aim of this work is to illustrate the use of the nonsymmetrized basis to deal 
with permutational-symmetry-breaking term in the Hamiltonian, in the present case the Coulomb 
interaction. 

I. INTRODUCTION 

The Hyperspherical Harmonics (HH) basis set has been extensively used to describe bound states and scattering 
processes in A = 3,4 systems One of the main reasons of using the HH basis resides in its flexibility, however 
it suffers from a large degeneracy, which, up to now, has prevented a systematic use with realistic potentials for 
A > 4 systems. The authors have recently proposed an approach 043 m which the HH basis is not symmetrized. 
In particular, in Ref. Q the nonsymmetrized basis has been used to describe systems up to A = 6 particles. It 
was shown that the large degeneracy of HH can be tackled noticing that the Hamiltonian can be expressed as an 
algebraic combination of sparse matrices, and the diagonalization procedure was implemented by means of an iterative 
diagonalization, where only the action of the Hamiltonian on a vector is required. In this work we continue the study 
on the A = 6 system, interacting via a two-body Volkov potential acting in s-wave, analyzing how the introduction 
of Coulomb interaction breaks the original permutation symmetry Sq, in the case of two or three protons. 

II. THE METHOD 

Following Ref. @, we introduce the Jacobi coordinates, x 1; . . . , x^r, with N = A — l, and we "adapt" the coordinates 
to the particle pair defining x^ = — rj, where and r 3 are the Cartesian coordinates of particles i and j. 
From the Jacobi coordinates, we introduce the hyperspherical coordinates, that means an hyper-radius p and 37V — 1 
hyperangles = (x x , . . . , x N , (f) 2 , ■ ■ ■ , <Pn)- 



The HH functions with fixed angular momentum LM, y^(ft 1 ^), are defined as the eigenfunctions of the grand 
angular momentum. They are labelled by the grand angular quantum number K and a set of 3iV — 1 quantum 
numbers indicated by [K]. In Ref. [4| it was shown that the potential energy, at fixed values of the hyper-radius, can 
be written as 

ij ij 

where the potential matrix [Vi2(p)][_R"'][if] = (y^(^jv)l^(l, 2)\y^ ! (Qj^)) is a sparse matrix and the matrix Bfj M is 

given by a product of sparse matrices corresponding to the transformation of the HH vector defined on to the one 
defined on The potential energy matrix is obtained after integrating on p using a Laguerre basis. Furthermore, 
the kinetic energy is diagonal in HH basis, and this allows to write the full Hamiltonian as an algebraic combination 
of sparse matrices. The diagonalization of the Hamiltonian is obtained via an iterative scheme, where only the action 
of the Hamiltonian on a vector is required. 
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TABLE I: Binding energies calculated with Volkov's potential in s-wave, for A = 6 particles and L = 0, with and without 
Coulomb interation, using 11 Laguerre's polynomials. 
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III. RESULTS 

For our numerical application, we chose a Volkov potential 

V(r) = V R e~ r2 / R * + V A ^ jRl , (2) 

with Vr = 144.86 MeV, Ri = 0.82 fm, V A = -83.34 MeV, and R 2 = 1.6 fin, which only acts in s-wave, and with the 
mass such that h 2 /m = 41.47 MeVfm 2 . We have calculated the binding energy for a system of A = 6 particles, and 
we repeated our calculations for the same system with Coulomb interaction between two particles, i.e. a model of 
6 He, and between three particles, i.e. 6 Li. 

In TableUwe show the results for the L = state. Without the Coulomb interaction the symmetry group is 56 and 
to antisymmetrize the wave function, taking also into account the spin and isospin degree of freedoms, the eigenvalue 
must belong to the irreducible representation [4 2]. We want to stress the fact that, using the nonsymmetrized basis, 
the eigenvectors belong to all of the irreducible representations of <S6, and in particular that of interest. 

When we add Coulomb interaction between two particles, the symmetry is broken as Sq — > S2 ® 5*4, and the original 
level, having degeneracy 9, is split into 4 sub- levels 

[4 2] -»• [l 2 ] <g> [3 1] + [2] ® [4] + [2] ® [3 1] + [2] <8 [2 2 ] , (3) 

where only [2] ® [2 2 ] is physical, as is the only one that can be antisymmetrized with respect the four neutrons using 
the spin degree of freedom, and describes the ground state of 6 He. When the Coulomb interaction is extended to 
three particles, the symmetry breaking is Se —> S3 ® S3, and the split reads 

[4 2] -+[21]® [3] + [21]® [2 1] + [3] <g> [3] + [3] ® [2 1] , (4) 

with the only physical state, describing 6 Li, being [2 1] ® [21]. 

We have shown the power and the flexibility of the nonsymmetrized HH approach. In particular we were able to 
include basis states up to K = 22 for a six-body system (corresponding to a basis set of 38 798 760 elements using 
11 Laguerre polynomials ). Furthermore, using a symmetry- adapted Lanczos algorithm, we were able to trace the 
irreducible representation of the eigenvector of interest and to select the corresponding eigenvalue. 
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